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ABSTRACT 

We present the first results obtained using a tunable narrowband filter in the search for high-z 
protoclusters. Using the recently commissioned red tunable filter on the Gran Telescopio Ca- 
narias we have searched for Lya emitters in a 75 arcmin 2 field centered on the z = 4.413 radio 
galaxy 6C0140+326. With three different wavelength tunings we find a total of 27 unique can- 
didate Lya emitters. The availability of three different wavelength tunings allows us to make 
estimates of the redshifts for each of the objects. It also allows us to separate a possible pro- 
tocluster from structure in the immediate foreground. This division shows that the foreground 
region contains significantly fewer Lya emitters. Also, the spatial distribution of the objects 
in the protocluster field deviates from a random distribution at the 2.5cr level. The observed 
redshift distribution of the emitters is different from the expected distribution of a blank field 
at the ~ 3cr level, with the Lya emitters concentrated near the radio galaxy at z > 4.38. The 
6C0 140+326 field is denser by a factor of 9 + 5 than a blank field, and the number density 
of Lya emitters close to the radio galaxy is similar to that of the z ~ 4.1 protocluster around 
TN J1338-1942. We thus conclude that there is an overdensity of Lya emitters around the 
radio galaxy 6C0140+326. This is one of few known overdensities at such a high redshift. 

Key words: galaxies: evolution - galaxies: high-redshift - galaxies: clusters: individual - 
cosmology: observations - cosmology: early Universe 



1 INTRODUCTION 

The identification of the progenitors of local galaxy clusters at 
z > 2 is a difficult task. For the interval 1 < z < 1.5, galaxy 
clusters are most often identified by infrared red sequence searches 
or observations of the X-ray emitting intracluster gas and the 
num ber of galaxy cluster s at these redshifts is growing steadily 
(e.g. IStanford et alj [l997l: iRosati et al.1 1 19991. 12004 iMullis etal 



20051: IStanford et all 120061 : iMuzzin et alj 120091 ; iBrodwin etal. 
20ialBielbvetalJl2010l) . Unfortunately, these methods become 
increasingly less effective when moving beyond z = 1.5 as the 
number of red galaxies decreases and X-ray emission becomes 
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too faint to be easily observed. However, if we wish to under- 
stand the role of environment on galaxy evolution and the emer- 
gence of large scale structure it is essential to locate and study 
galaxy clusters at all possible epochs. Recent results have pre- 



ray emission at z > 1.5 dWilson et alj 


20081; iKurk et al. 


2009; 


Paoovich et al. 2010: Tanaka et al. 20ld: 


Henrv et alj 201C 


), with 



the current distance record bei ng the galaxy c luster CLJ 1449+0856 
at z ~ 2.07 presented by iGobat et alj d201 lh (but see also 
lAndreon & Huertas-Companvl 1201 If) . However, this sample of 
high-z clusters remains small. 

Another successful method of identifying galaxy cluster 
progenitors at z > 2 is to search for overdensities of line 
emitting galaxies using narrowband imaging. These searches 
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are often aimed at fields containing high-z radio galaxies 
(hereafter HzRGs, iMilev & De Breuc kl l2008t) . since these are 
thought to have large stellar masses of the order of 10 " to 
10 12 M Q jRocca-Volmerange et al.l2004ISevmour et al.l2007l) . Ac- 
cording to hierarchical galaxy formation, the most massive galaxies 
form in the densest environments. The massive nature of HzRGs 
thus indicates that these objects may trace overdensities in the 
early Universe. In recent years many studies h ave focused on 



findin g galaxy overdensiti es aro und HzRGs (e.g . Pascarell e et al 



199d Rnopp & Chamberslll997l;|Pentericci et : alj|2000l;lKurk etal 



2004a. b; Overzier et al. 2006; Venemans et al. 2007; Overzier et al 
20081; Kuiper et alJl2010l : lGalametz et al]l2010r . lHatch et alfeouf) . 
Since these overdensities show no evidence of X-ray emission (at 
luminosities > 10 44 erg s~') it is t hought that these are forming clus- 
ters, n ot yet dynamically relaxed 1 Car illi et al.l2002l ; IOverzier et al.l 
2005). They are therefore often called 'protoclusters'. 

However, even though the number of spectroscopically con- 
firmed z > 2 HzRGs approaches 200, the fraction of these that 
have been studied for the presence of galaxy overdensities remains 
small. This is partly due to the small number of existing narrow- 
band filters and the fact that the central wavelengths of these ex- 
isting filters are often based on strong lines at z = 0, such as 
[Om]/15007. This severely limits the redshifts at which protoclus- 
ters can be studied and therefore the absolute number of confirmed 
protoclusters has remained small. 

In this paper we present the results of a pilot study that utilizes 
tunable narrowband filters in the search for line emitting galaxies 
around the HzRG 6C0 140+326 at z = 4.413 dRawlings et all 19961 ; 
|Pe Breuck e"tafll200ll) . Tunable filters (TFs) allow the user to set 
the central wavelength and width of the narrowband filter. TFs use 
two plane parallel transparent plates coated with films of high re- 
flectivity and low absorption. By separating the two plates by a 
small distance of the order of a yum-mm a cavity is formed which 
is resonant at a specific wavelength. Constructive interference at 
the resonant wavelength then causes all the incident light at that 
wavelength to be transmitted. Changing the separation between the 
plates then allows the central wavelength or the width of the filter 
to be adjusted. TFs therefore alleviate the limitations imposed by a 
small number of available narrowband filters at fixed wavelengths 
and are thus ideally suited for searching for protoclusters at a range 
of redshifts. Mor e information and detail s co ncerning tunable fil ters 
can be found in lBland-Hawfhornl dl995h and ljones etal] J2002h . 

Similar studies involving the search for line emitting galax- 
ies around z ~ 1 quas ars have been su ccessfully performed by 
iBaker et all fcOOlh and iBarr et"ai] d2004h . An attempt to use this 
technique at hig her redshifts has led to mixed results. In unpub- 
lished work by ISwinbankl fc006t) a TF study is presented for two 
radio-loud quasars located at z ~ 2 and one radio-quite quasar 
at z ~ 4.5. The z ~ 4.5 field shows evidence for an overdensity, 
whereas the two z ~ 2 fields lack depth and do not allow for strong 
conclusions. With the advent of a TF instrument at a 8-10 meter 
class telescope, it has now become possible to obtain sufficiently 
deep data to efficiently search for Lya emitters in the environments 
of HzRGs at arbitrary redshifts z > 2. 

The paper is organised as follows: in Sect. [2] we describe the 
data, its reduction and the object detection. The sample selection 
and redshift estimation is treated in Sect. [3] and we discuss the ev- 
idence for the presence of an overdensity in Sect. [4] Finally, con- 
clusions and future outlook are presented in Sect. \5\ Throughout 
this paper we use a standard ACDM cosmology with H = 71 km 
s Mpc~', Q. M = 0.27 and S1 A = 0.73. All magnitudes given in this 
paper are in the AB magnitude system. 



2 DATA 

The radio galaxy 6C0140+326 (hereafter 6C0140) was observed 
for a total of 18 hours using the Optical System for Imaging 
and low Resolution Int egrated Spectroscopy instrument (OSIRIS, 
ICepa et alj[2000L 120031) at the Gran Telescopio Canarias (GTC), 
La Palma. OSIRIS consist of two 2048 x 4096 pixel Marconi 
CCDs with a 72 pixel gap between the two CCDs. The observa- 
tions were done on several dates from September 2010 to January 
2011. The 2x2 binning mode was used resulting in a pixel scale 
of ~ 0.25" pixeL 1 and a total field-of-view of ~8.7'x8.6'. The ra- 
dio galaxy was positioned near the optical centre, approximately 
15" from the left edge of CCD 2. The individual exposures have 
been dithered with offsets of approximately 10-12.5" in right as- 
cension and 2-4" in declination. The larger offsets in right ascen- 
sion have been chosen such as to cover the gap of ~ 8 " between 
the two CCDs without losing the radio galaxy in the gap. A full list 
of details concerning the observations can be found in Table Q] 

Broadband images were obtained in the r and ( bands. The 
narrowband images were obtained at three different central wave- 
lengths /l c . This was done because the central wavelength of the 
TFs varies across the field of view approximately as 



A(r) = A c (l - 0.0007930/" 2 ) 



(1) 



as given in the OSIRIS TF user manual. Here r is the distance to 
the optical centre of the instrument in arcminutes and A c the wave- 
length at the optical centre. The maximum FWHM of the TF is 
20 A, therefore it was necessary to perform several passes with dif- 
ferent A c in order to cover the entire redshift range of a possible pro- 
tocluster. At z ~ 4.4 the Lya- line is shifted to 6580 A, so the central 
wavelengths were chosen to be 6565 A, 6575 A and 6585 A (here- 
after TF1, TF2 and TF3 for brevity). Note that the redshift of the 
radio galaxy (z = 4.413) indicates that it falls between TF2 and 
TF3. 

As noted above, the maximum formal width of the TF is 20 A. 
However, the shape of the response curve is Lorentzian rather than 
Gaussian and is approximately given by 



1 + 



2(A - A c ) 



6A 



(2) 



with A c the central wavelength and SA the formal FWHM. Due 
to the shape being Lorentzian, the transmission has relatively ex- 
tended wings which results in an effective band width that is larger 
than the formal value for the FWHM by a factor of n/2. Thus the 
TF tunings have an effective FWHM of ~ 31 A. Taking this into ac- 
count, our observations probe the redshift range 4.386 < z < 4.428 
near the optical centre. The true redshift range that is covered is 
larger due to the variation of the central wavelength across the field. 
The relevant filter response curves are shown in Fig.Q]together with 
a night-sky emission line spectrum. As can be seen, the sky line 
contamination is relatively mild. 



2.1 Data reduction 

The reduction of all the data is done using ira^. The reduction 
of the broadband images includes the standard steps of bias sub- 
traction and flat-fielding, where the latter is done using sky flats. 



1 iraf is distributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 
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Table 1. Details of the observations. The 5<x limiting magnitudes have been calculated for an aperture diameter of twice the seeing disk. Also note that the 
values given for A^g for the TF observations are the values at the optical centre. 



Band 


Exp. time (sec.) 


ieff (A) 


AA (A) 


Seeing (arcsec) 


5<r limiting magnitude 


r 


2400 


6417 


1685 


0.8 


25.8 


i 


2100 


7719 


1483 


0.8 


24.8 


TF1 


15680 


6565 


31.4 


0.8 


25.1 


TF2 


15680 


6575 


31.4 


0.8 


25.2 


TF3 


13440 


6585 


31.4 


0.8 


25.1 




5500 6000 6500 7000 6350 6400 6450 6500 6550 6600 6650 6700 

Wavelength (angstrom) Wavelength (angstrom) 

Figure 1. Left panel: Filter response curves for the three TF tunings (solid curves) and the r band (dashed curve). Right panel: A zoom in of the three TF 
response curves in order to show more detail. This panel shows also the approxima te area wa velength range where protocluster galaxies are expected to lie. 
Also shown in grey in both panels is the night-sky emission line spectrum (Hanuschik 2003). The minimum wavelength reached at the edge of the field in 
TF1 is denoted by the vertical dotted line. 



To remove further large scale gradients a superflat is made from 
the unregistered science images. This superflat is smoothed and the 
science frames are subsequently divided by the superflat. The flat 
science images are then registered with simple offsets in the x and 
y direction using the iraf task xregister and co-added together. 

The reduction of the TF images follows the same general out- 
line as that for the broadband images with two exceptions. The flat- 
fielding is done with dome flats and an additional step is included 
which involves the removal of sky rings. As the central wavelength 
of the TF filter varies across the field, skylines shift in and out of 
the filter bandpass causing a pattern of alternating bright and faint 
concentric rings superimposed on the image. These large scale gra- 
dients are of the order 3 to 9 times the rms noise in the individual 
images. The rings are removed by subtracting a smoothed superflat 
made using the unregistered science exposures. This superflat is 
created for each exposure individually as the sky level (and there- 
fore the sky rings) varies between different nights and airmasses. 
The individual frames are subsequently registered and combined. 
Finally, all fully reduced science images are registered to the same 
pixel coordinates using the iraf tasks geomap and georegister. 

Due to the wavelength variation across the field and the dither- 
ing there is a variation in wavelength in each of the individual 
pixels. This variation is larger near the edge of the field of view. 
The wavelength assigned to each pixel in the final images is the 
mean of the wavelengths of the pixel in question in the individ- 
ual images. This also implies that the effective FWHM of the TF 
increases when moving away from the optical centre. This effect 
is strongest in the right ascension direction because the dither- 



ing steps are larger in this direction, with a maximum increase of 
~ 50 per cent at the very edges of the images. 

Flux calibration for both broad and narrowband imaging is 
achieved using standard star observations obtained at the end of 
each observing block of 1 hour. The standard stars used for the TF 
flux calibration have full SEDs available allowing flux calibration 
at the exact wavelength of each of the TF observations. A further 
independent check of the flux calibration of the TF observations is 
obtained using the broadband data. The r and i magnitudes of all 
objects in the science frames with 19 < r < 23 are measured. Then, 
assuming a power law spectral energy distribution, the magnitudes 
of these objects at the wavelength of interest are determined. The 
median zeropoints derived with this method deviate by -0.1, -0.04 
and +0.005 magnitude with respect to the standard star zeropoints 
of TF1, TF2 and TF3. The larger deviation for TF1 and TF2 are 
likely due to the stronger presence of the Ha absorption line at 
6563 A. 

The final reduced and coadded TF images show some arti- 
facts of the reduction. The sky ring subtraction is not optimal due 
to the applied smoothing and therefore a residual ring pattern re- 
mains in the final images. Also the unique properties of the TF lead 
to pupil ghosts near bright stars. Point source ghosts, however, are 
not present in the final TF images. The dithering results in an offset 
in point source ghosts opposite to the actual offset. When combin- 
ing the individual images the point source ghosts will therefore be 
removed. 
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2.2 Source detection and photometry 

Object detection and photometry are done using SExtractor 
jBertin & Arnoutslll996h in double image mode. We create three 
different catalogues based on each of the three TF tunings. A de- 
tection is defined as a minimum of 9 adjacent pixels that each 
exceed the 2<r rms noise. Colours are measured using the 2cr 
isophotal apertures as determined from the respective detection im- 
age whereas total magnitudes are measured using SExtractor's 
mag_auto apertures. Image depth and u ncertainties on the p hotom- 
etry are determined using the method of lLabbe et al.l (2003). 

Completeness of the image is measured by adding point 
sources of a range of magnitudes to the respective images after 
which source extraction is repeated and the number of recovered 
objects is assessed. To avoid overcrowding the image we limit the 
number of added objects to 150 per magnitude. This process is re- 
peated 10 times in order to obtain better statistics. We find that 
the data is 50 per cent complete for point sources down to r=26.0, 
TF1=25.1, TF2=25.1 and TF3=25.0 magnitudes. 



3 RESULTS 

3.1 Selection of LAEs 

For each of the three TF tunings a separate sample of Lya emitters 
(hereafter LAEs) is selected. The criterium for identifying LAEs 
is based on the colour-magnitude diagrams of the 6C0140 field as 
shown in Fig. [2] The galaxies with line emission at the relevant 
redshift will show an excess of flux in the TF band relative to the 
r band flux. Most objects do not have emission lines in the TF and 
thus have r-TF ~ 0. The rms scatter around r-TF ~ at TF ~ 24- 
25 mag is ~ 0. 15. For an object to be identified as LAE we require 
at least a 5cr deviation from r - TF ~ and thus that r - TF > 0.75. 

This excess flux relates to an observed line equivalent width 
using the relationship 

AVLItp[i - 10-° 4 <'- TF >] 

EWob * = [AA r W-o^-^ - AA TF ] (3) 

from lBunkeretal.1 jl995h . Here AA is the FWHM of the filter in 
question and r and TF are the measured magnitudes in the broad- 
and narrowband respectively. The restframe equivalent width EWq 
is obtained by dividing by (1 + z). Note that this relation does not 
account for IGM absorption. For LAEs at z ~ 4.4 this implies that 
EW \, S and EWq are overestimated by a factor of ~ 1.5-2. We will 
return to this correction factor in Sect. |3.2| 

Using Eq. [3] the colour cut corresponds to approximately 
EW \, S > 32 A at the optical centre of the field. However, as de- 
scribed in Sect. 12. 11 the dithering results in an effective broadening 
of the TF when moving away from the optical centre. At the West 
and East edges of the image the increase in effective FWHM is 
~ 50 per cent, thus the equivalent width cut varies approximately 
between EW obs > 32 A and EW obs > 48 A across the field. 

The equivalent width cut used in this work is less stringent 
than the more commonly used EW obs > 80 A. This is because the 
TF is narrower than conventional narrowband filters, thus allowing 
for lower equivalent width objects to be included. We further re- 
quire that the error parameter X > 3, i.e. the excess flux is at least 
three times larger than the combined noise of the measured broad- 
and narrowband fluxes. Finally, due to the artifacts present in the 
images we visually inspect all of the objects that satisfy the above 
criteria and discard any spurious detections. 

As can be seen from Fig. [2] there is a significant number of 



objects located in the selection area that are not identified as be- 
ing LAEs. Visual inspection shows that these objects are likely 
spurious detections because they are often found in the sky ring 
residuals, have unphysical shapes, coincide with the pupil ghosts 
or are located near bright stars and saturation spikes. These objects 
are therefore not included in the candidate LAE sample. Figure [3] 
shows the regions of the image where most of these spurious are 
located, i.e. the sky rings and the pupil ghosts. Also shown are the 
locations of bright stars and portions of the image that are affected 
by vignetting. The bright stars cover ~ 4 per cent of the field and 
therefore do not influence any of the conclusions presented in this 
work. 

Other objects close in the selection area may actually have 
E < 3. Individual values of E depend on the aperture size for each 
object. The 2 = 3 curves shown in Fig. [2] however, are calculated 
for a fixed aperture size that is taken to be the median aperture size 
of the LAEs. Thus individual objects inside the selection area may 
have X < 3 or vice versa. 

For the TF1 filter, which probes exclusively z < 4.407, we find 
a total of 13 candidate LAEs, including the radio galaxy which has 
EW obs = 88.3 A. This is smaller than expected as visual inspection 
of the TF1 image indicates that the radio galaxy does have a strong 
narrowband excess. However, the Lya emi ssion is extended an d 
coincides with a z ~ 0.9 foreground galaxy (Rawlinas et alJI 19961) . 
This results in strong contamination of the r band flux and therefore 
the expected narrowband excess is significantly diminished. To al- 
leviate this problem we use the r band to define the colour apertures 
for the radio galaxy. This results in EW bs = 217 A indicating that 
the presence of the foreground galaxy is indeed important. 

The TF2 band at 6575 A reveals a total of 14 candidate LAEs, 
including the radio galaxy with EW obi = 470 A. Eight of the LAEs 
have been identified in TF1 as well. 

Finally, a total of 17 objects are identified as being candidate 
LAEs in the TF3 band at 6585 A. Again, the radio galaxy has been 
included and its equivalent width is highest in TF3 with EW obs = 
578 A. Out of the remaining galaxies eight have not been detected 
in TF1 or TF2, thus yielding a total number of 27 unique LAEs 
(including the radio galaxy) in the field of 6C0140. The properties 
of the emitters are listed in Tableland Table[3] 



3.2 Redshift distribution 

Having the multiple tunings of the TF means the redshift of the 
LAEs can be constrained to a greater accuracy than with a single 
narrowband image. This allows us to investigate the approximate 
redshift distribution of the candidate LAEs in the field. 

In order to obtain the best possible estimated redshifts we 
model what the effect is of a variety of weighting schemes. We 
model a simple flat spectrum (/3 = -2) with an emission line with 
restframe FWHM of 250 km s _I or 500 km s -1 at certain wave- 
lengths in the range 6555 < A < 6595 A. These values for the 
FWHM are consistent with the values found for 80 to 90 per cent of 
LAEs around HzRGs ( tVenemans et al. 2005 1120071). A bsorption by 
the IGM is taken into account using the lMadaul dl995h recipe. The 
modeled spectra are then convolved with the TF response curve as 
given in Eq.[2]with A c = 6565, 6575 and 6585 A. Using the fluxes in 
each of the TF tunings we then derive estimates for the wavelength 
of the emission line. 

The results of this process are shown in Fig. [4] We show two 
different weighting schemes which can be given in general form as 




Figure 2. Colour-magnitude diagrams for each of the three TF tunings. The horizontal solid lines indicate EW b s values of 32 and 80 A respectively, whereas 
the dashed curve indicates the S = 3 line. Objects identified as LAEs are denoted by diamonds and labeled with the respective ID numbers as listed in Table|2] 
and Table[3] The radio galaxy is denoted by the square symbol. The X = 3 curve is calculated using a median aperture size. Individual objects may therefore 
be located inside the selection area as depicted here, but are still not identified as being a LAE. These objects are denoted by open circles. Spurious detections 
are indicated by crosses. 



Table 2. Coordinates and magnitudes of the detected LAEs. The table is divided in three sections. The objects in the first section are brightest in TF1 and so 
forth. "The radio galaxy. The r band was used as the detection image for this object. 



ID 




RA 


Dec. 


Detected in 


r 


TF1 


TF2 


TF3 


#1 


01 


43:58.2 


+32:49:49.8 


TF1 


27.2 ±0.8 


24.9 ± 0.2 


25.34 ± 0.3 


26.3 ±0.6 


#2 


01 


43:44.6 


+32:49:42.6 


TF1 


23.1 ±0.1 


21.87 ±0.04 


22.55 ± 0.06 


22.9 ±0.1 


#3 


01 


43:44.8 


+32:57:06.0 


TF1,TF2,TF3 


26.7 ± 0.4 


24.7 ±0.1 


25.0 ± 0.2 


25.4 ± 0.2 


#4 


01 


43:38.9 


+32:53:37.2 


TF1,TF2 


> 27.5 


24.9 ± 0.2 


25.2 ±0.2 


25.7 ±0.3 


#5 


01 


43:38.3 


+32:49:58.3 


TF1 


24.6 ±0.1 


23.7 ±0.1 


23.9 ±0.1 


24.2 ± 0.2 


#6 


01 


43:35.7 


+32:54:36.0 


TF1 


25.7 ±0.2 


24.7 ±0.1 


25.1 ±0.2 


25.4 ±0.3 


#7 


01 


43:35.4 


+32:52:11.5 


TF1 


26.7 ± 0.3 


25.3 ±0.2 


25.8 ±0.2 


25.8 ± 0.2 


#8 


01 


43:24.8 


+32:54:09.5 


TF1/TF2 


26.5 ± 0.4 


24.3 ±0.1 


25.4 ±0.3 


25.1 ±0.2 


#9 


01 


43:45.9 


+32:53:58.0 


TF1.TF2 


25.9 ± 0.2 


24.4 ±0.1 


24.4 ±0.1 


24.8 ± 0.2 


#10 


01 


43:43.3 


+32:52:08.1 


TF2.TF3 


27.4 ± 1.0 


25.5 ± 0.3 


24.6 ±0.1 


24.8 ± 0.2 


#11 


01 


43:42.1 


+32:55:27.6 


TF2 


23.7 ±0.1 


23.1 ±0.1 


22.8 ±0.1 


23.1 ±0.1 


#12 


01 


43:41.5 


+32:54:17.2 


TF1,TF2,TF3 


26.6 ±0.5 


25.3 ± 0.3 


24.3 ±0.1 


24.5 ± 0.2 


#13 


01 


43:41.4 


+32:53:49.3 


TF1,TF2,TF3 


25.6 ±0.3 


24.6 ± 0.2 


24.0 ±0.1 


24.4 ±0.1 


#14 


01 


43:41.0 


+32:53:48.5 


TF2 


27.2 ± 0.5 


26.0 ± 0.3 


25.7 ±0.2 


26.3 ± 0.4 


#15 


01 


43:40.2 


+32:55:02.5 


TF2.TF3 


25.9 ±0.4 


24.7 ± 0.2 


23.9 ±0.1 


24.1 ±0.2 


#16 


01 


43:38.0 


+32:49:51.9 


TF1,TF2,TF3 


22.9 ±0.1 


22.19 ±0.07 


21.85 ±0.05 


22.35 ± 0.08 


#17 


01 


43:59.8 


+32:52:16.2 


TF3 


25.1 ±0.1 


25.0 ± 0.2 


24.8 ±0.1 


24.2 ±0.1 


#18 


01 


43:56.4 


+32:54:41.1 


TF2.TF3 


> 27.5 


26.2 ± 0.3 


25.7 ±0.2 


25.5 ±0.2 


#19 


01 


43:44.8 


+32:56:01.9 


TF3 


26.7 ± 0.2 


26.7 ± 0.4 


27.0 ± 0.5 


25.3 ±0.1 


#20 


01 


43:43.9 


+32:52:28.6 


TF3 


> 27.5 


27.5 ± 1.6 


25.7 ±0.3 


24.9 ±0.1 


#21" 


01 


43:43.8 


+32:53:49.9 


TF1,TF2,TF3 


23.65 ± 0.05 


21.54 ±0.01 


20.91 ± 0.01 


20.75 ± 0.01 


#22 


01 


43:38.0 


+32:52:00.9 


TF3 


25.8 ±0.1 


25.8 ± 0.2 


25.6 ± 0.2 


24.9 ±0.1 


#23 


01 


43:36.1 


+32:55:00.9 


TF2.TF3 


24.0 ±0.1 


23.8 ±0.1 


23.17 ±0.05 


23.13 ±0.05 


#24 


01 


43:33.3 


+32:54.09.0 


TF3 


25.6 ±0.2 


26.1 ±0.5 


25.2 ± 0.2 


24.6 ±0.1 


#25 


01 


43:27.3 


+32:51:32.4 


TF3 


27.1 ±0.5 


27.5 ± 1.4 


25.7 ±0.3 


25.1 ±0.2 


#26 


01 


43:26.6 


+32:52:17.3 


TF3 


25.7 ± 0.2 


27.1 ±0.9 


25.5 ±0.2 


24.8 ±0.1 


#27 


01 


43:25.9 


+32:52:01.8 


TF3 


26.0 ± 0.2 


26.4 ± 0.6 


25.9 ±0.3 


25.1 ±0.2 



Here i is TF1, TF2 or TF3 respectively, F t is the flux in the respec- 
tive bands, A t is the central wavelength of the tuning in question at 
the relevant location in the field and w, is a weighting factor. The 
black data points represent the case of w, = 1 for all bands and the 
red data points denote the results obtained for 

w . — 2 f ';/ rain ( f TFl.fTF2.' r TF3) (g\ 



As noted above, the results shown in Fig.[4]are obtained when us- 
ing the central wavelength values of the optical centre, i.e. Ajfi = 
6565 A, Ajp2 = 6575 A and /i TF3 = 6585 A. However, the wave- 
length shift acts on each tuning identically and thus the qualitative 
behaviour shown in Fig. [4] is valid for the entire field, irrespective 
of wavelength shift. 

We see that the more involved weighting scheme yields bet- 
ter results for almost all redshifts. However, for both weighting 
schemes the largest discrepancy between the input and output red- 
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Figure 3. A mask of the 6C0140 field showing the location of regions with 
large numbers of spurious detections and regions where source detection is 
impeded. The large grey and black rings are the approximate locations of 
the sky rings in TF1 and TF3, respectively. For TF2 the sky rings fall in 
between the two sets shown here. The small black arcs denote the locations 
of the pupil ghosts, whereas the filled circles denote the location of bright 
stars. The two light grey regions at the left side of the image are affected by 
vignetting and therefore not considered. 



shifts is at either end of the investigated redshift range. This is to be 
expected, because no data is available to bracket the existing tun- 
ings, hence skewing the output redshift towards a central value. The 
more elaborate weighting scheme alleviates this slightly, but does 
not yield full agreement. Figure|4]also shows that the results of the 
two weightings schemes are fairly robust with respect to the choice 
of restframe FWHM. 

The redshifts obtained when using the weights as described 
in Eq. [5] are listed in Table [3] Uncertainties have been calculated 
by varying the measured fluxes according to their respective un- 
certainties and recalculating A,.g- An additional systematic uncer- 
tainty was added in quadrature to take into account the limitations 
of the weighting scheme. This systematic uncertainty is based on 
the bands in which the object is detected and the offset from the 
input redshift implied by this as measured from Fig. [4] Finally, to 
account for the variation in wavelength in each pixel we add an 
additional location-dependent uncertainty in quadrature. 

None of the LAEs are located at larger redshifts than that of 
the radio galaxy. This is a selection effect as the highest value of A c 
is 6585 A whereas the Lya line for z RG =4.413 falls at ~ 6580 A. 
Combined with the wavelength shift towards shorter wavelengths 
across the field this means that the observations are biased towards 
redshifts lower than the redshift of the radio galaxy. Also, note that 
the estimated redshift of the radio galaxy is only marginally in- 
consistent with the spectroscopic redshift of z = 4.413. The slight 
underestimate with respect to the spectroscopic redshift is possi- 
bly due to the fact that we modeled the Lya line with a maxi- 
mum FWHM of 500 km s~' , whereas the line is observed to have a 
FWHM of ~ 1500 km s~' (Rawlings etal.ll 19961 : 1 De Breuck et al. 
l200ll) . A larger line width will introduce stronger systematic un- 




4.395 4.400 4.405 4.410 4.415 4.420 



Figure 4. Upper panel: Output versus input redshift when using the two dif- 
ferent weighting schemes tested on modeled Lya spectral lines. The solid 
curves indicate the results obtained when using = 1, whereas the dia- 
monds denote the results obtained with the weights described as in Eq.[5] 
Black and red indicate restframe FWHM of 250 km s -1 and 500 km s _I , re- 
spectively. Lower panel: Difference between the output and input redshifts. 

Table 3. Properties of the detected LAEs. "EW \> S calculated usin g Eq. [3] 
whereas EWq is calculated using the method of Venemans et al. (2005). 
b EW a \ >s , EWo, £ and F^ ya for the radio galaxy are calculated on the ba- 
sis of TF3. The horizontal dividers are as in Tablef2] 



ID EWobs/EWo (A)° 2 F Lya (erg s" 1 cnT 2 ) 



#1 


411/36.2 


5.2 


1.4 x 10~ 


17 


4 o| +0.005 


#2 


84.5/9.1 


17.9 


9.0 x 10~ 


17 


4 ^7+0.002 

-Q-QP.S 


#3 


209/24.7 


6.8 


1.7 x 10~ 


17 


4 3 <r 9 +8.883 
-Q.QQ7 
4 3Q4+O.O85 


#4 


> 463/ > 49.8 


6.2 


1.7 x 10- 


17 


#5 


47.6/8.0 


5.0 


2.3 x 10- 


17 


4 340+ - 004 
-P..QA9 


#6 


60.3/6.5 


4.8 


1.5 x 10- 


17 


4380 -: 


#7 


96.9/13.0 


4.7 


1.8 x 10- 


17 


4 377+8.083 

' -0.009 


#8 


266/29.6 


8.9 


2.2 x 10- 


17 


4 "577+0.002 
-0.008 


#9 


106/11.7 


6.7 


1.9 x 10- 


17 


4 4 n 4 +0.002 


#10 


552/61.5 


6.9 


1.5 x 10- 


17 


4 399+0.004 

4 391+8:881 

*- Jy -0.002 


#11 


47.5/5.0 


10.3 


3.7 x 10- 


17 


#12 


299/31.0 


7.4 


2.5 x 10- 


17 


4 404+0- 004 


#13 


122/12.9 


7.2 


3.7 x 10- 


17 


4 404+0-002 


#14 


120/12.4 


3.5 


7.0 x 10- 


18 


4 403+ - 003 


#15 


225/23.5 


7.8 


4.3 x 10- 


17 


4395+0.003 

4 340+8:883 
* -0.003 


#16 


68.8/7.1 


12.3 


1.1 x 10- 


16 


#17 


52.0/4.5 


5.8 


1.1 x 10- 


17 


4 ^ftR+O- 009 


#18 


> 247/ > 25.9 


5.0 


1.1 x 10- 


17 


4 +8.088 
-0.006 


#19 


96.8/10.0 


6.0 


5.5 x 10- 


18 


4 QO7+8.009 
t.jo;_ 004 


#20 


> 434/ > 45.0 


6.4 


1.9 x 10- 


17 


4 40q+ - 008 


#21* 


578/75.0 


95.5 


3.6 x 10- 


15 


4419+O.OOI 
H -^ lz -0.001 


#22 


47.3/4.5 


5.0 


6.9 x 10- 


18 


4 394+O.OO8 
-Q-QQ3 


#23 


45.6/4.6 


8.7 


2.8 x 10- 


17 


4-388™ 


#24 


65.5/5.9 


4.8 


1.6 x 10- 


17 


4-388!™ 


#25 


300/24.3 


5.2 


1.8 x 10- 


17 


4 339+0.009 


#26 


68.2/5.8 


5.3 


9.6 x 10- 


18 


4 344+8:885 


#27 


57.8/5.3 


3.2 


1.0 x 10- 


17 


4 ^q+Sffi 
■ -0.005 
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certainties that have not been taken into account in the case of the 
estimated redshift of the radio galaxy. This may therefore account 
for the discrepancy. 

The redshift estimates can be used to correct the equivalent 
width values for IGM absorption. As discussed in Sect. 13. II Eq.[3] 
does not take this into account and underestimates the r contin- 
uum flux density. Therefore EW h* and EWn are overestimated. The 

I ; 1 

correc tions are calculated following the method of Venemans et al. 
(2005). The resulting corrected restframe equivalent widths and 
corresponding Lyo- fluxes are listed in Table [3] The difference be- 
tween EWobs and EWo is typically a factor ~ 10. This is consistent 
with a factor (1 + z) in combination with a factor ~ 1.7. Here the 
latter factor originates from the fact that a larger portion of the r 
band flux is absorbed by the IGM compared to the TF flux. The ex- 
act factor varies between ~ 1 - 2 and depends on the location of the 
line within the filters and therefore both the redshift of the object 
and its position in the field. 



3.3 Contamination 

One of the larger caveats of using narrowband imaging to select 
high-z emission line galaxies is that the final sample may be con- 
taminated by low-z interlopers that have a strong emission line 
falling in the narrowband. For our study the most likely interlopers 
are [On] emitters at z ~ 0.76 or [Om] emitters at z ~ 0.3. Spectro- 
scopic follow-up is needed to accurately determine the success rate 
of the sample presented in this work. However, based on previous 
spectroscopic studies of z ~ 4 narrowband surveys the expected 
number of interlopers in our sample can be estimated and it can be 
determined whether the results presented here are robust when this 
is taken into account. 

We base our estimate of the success rate on the studies of 
HzRG TN J1338-1942 (hereafter 1338) bvlVenemans et al.l J2002l. 
|2007|) . the field study at z ~ 4.5 by iDawson et alj d2007l) and 
the study of the LAEs around z ~ 5.2 HzRG TN J0924-2201 by 
Venemans et al. I J2004T) . The success rates in each of these works 
are fairly similar to each other, ranging from ~ 75 per cent to 
~ 95 depending on whether non-detections are counted as non- 
confirmations. To investigate the 'worst case scenario' we use for 
our sample the minimum success rate of 75 per cent. 



4 DOES 6C0140+326 RESIDE IN A PROTOCLUSTER? 

Based on our sample of LAEs we determine whether there is an 
overdensity around 6C0140. Due to the wavelength shift across the 
field, part of the observed field can act as a control field. 

In Fig. [5] the spatial distribution of the LAEs is shown. Also 
shown are two concentric circles indicating the boundaries of two 
fields: the inner circle of ~ 12.5 arcmin 2 and the annulus of 
~ 16.3 arcmin 2 , respectively. The central field covers the redshifts 
closest to the radio galaxy and can therefore be considered a pos- 
sible protocluster field. The annulus outside this field delimited by 
the dashed circle probes lower redshifts and is considered to be 
field environment. To make the distinction between protocluster 
and foreground stronger we will only consider TF2 and TF3 de- 
tected objects (diamonds and asterisks) in the central field, whereas 
in the annulus only TF1 detected objects are considered (squares). 
This selection effectively means we are limited to z > 4.38 in the 
central field and z < 4.38 in the annulus. Furthermore, the width of 
the TF does not change significantly across the central field and the 
annulus. Since we only consider TF1 in the annulus, the physical 




□ 



□ 



25.95 

RA (degrees) 



Figure 5. The spatial distribution of the LAEs. Objects detected in TF1, TF2 
and TF3 are marked by squares, diamonds and asterisks respectively. Each 
object is also marked by its ID number. The location of the radio galaxy 
is marked by a star. Also shown are two circles denoting the border of a 
protocluster field (within the full circle) and a foreground field (between 
dashed and full circle). 



depth of the annulus is thus ~ 1.3 times smaller. We thus require the 
area of the annulus to be larger by the same factor to have the same 
volume in each of the fields. We find 9 objects (excluding the radio 
galaxy) in the possible protocluster field versus 2 in the foreground 
field. This thus indicates that there is a concentration of LAEs near 
to the radio galaxy. 

It is striking that almost all of the objects within the protoclus- 
ter field are located west of the radio galaxy in a North-South fil- 
amentary structure. We test whether the spatial distribution is con- 
sistent with a random distribution by applying a two dimensional 
Kolmogorov-Smirnov test. There is a probability of 0.01 that the 
distribution as shown in Fig. [5] is drawn from a random distribu- 
tion. The distribution is thus different from random at the ~ 2.5cr 
level. This further indicates that the LAEs are clustered. 

Since our control field is not very large, it is susceptible 
to cosmic variance. To better quantify the overdensity of LAEs 
around the radio galaxy, we also compare it to the blank field 



Dawson et 



also comp; 
afl j2007l) .1i 



Dawson etaljj2007l) pre- 



LAEs observed by 
sented a differential hya luminosity function for field LAEs with 
EW obs > 80 A at z ~ 4.5. Fitting a Schechter function to the 
luminosity function using a fixed value of a = -1.6 they find 
O* = (1.7±0.2)x 10~ 4 Mpc- 3 andL* = (10.9 ± 3.3) x 10 42 ergs" 1 . 

For the 6C0140 field, if we take into account the overlap be- 
tween the different tunings, the total 'unique' volume probed by the 
prot ocluster field is 1570 M pc 3 . Using the same selection procedure 
as in lDawson et aljj2007l) , we find eight emitters (excluding the ra- 
dio galaxy) in the central field. We calculate the expected number 
of LAEs in the same field (assuming it is a blank field) using the 
field luminosity function. The expected number of LAEs is found 
to be 0.9+j} 4 , with the uncertainty derived from the uncertainties 
on O* and L*. Here we use the Lya flux of the faintest emitter 
with EW obf , > 80 A (#14) in the protocluster field as lower limit 
(7.0 x 10~ 18 erg s -1 cirT 2 , this equals L = 1.4 X 10 42 erg s" 1 for 
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the cosmology used bv lDawson et alj2007l) . We therefore find that 
the 6C0140 field is denser than a blank field by a factor 9 ± 5, 
where the uncertainty is based on Poisson statistics and the un- 
certainty on the expected number. Defining galaxy overdensity as 
S g = «ciustcr/«fieid - 1 , we thus find a galaxy overdensity of 8 ± 5. 

Using the overdensity in 6C0140 protocluster field we can 
make a rough estimate of th e mass contained in this central field. 
As in lVenemans e t al. (2005) we use the relation 

M =P v [ l + jj (6) 

with p the mean density of the Universe, V the comoving volume 
considered and b the bias parameter which rel ates the galaxy over - 
den sity to the matter overd ensity. Following ISteidel et al.l d 19981) 
and Ishimasaku et al] {2003) we use b = 3 - 6. Using p = 3.5 x 
10'° M Q Mpc- 3 and V = 1570 Mpc 3 a mass of 0.8-2.9X 10 14 M Q is 
found. This is a strict lower limit to the mass of the entire overden- 
sity because the true extent of the protocluster is likely larger than 
what is indicated by the central protocluster field. 

We also determine whether the number density and redshift 
distribution of LAEs in the 6C0140 field are consistent with that of 
a z ~ 4 protocluster. However, the number of known protoclusters 
above z = 4 is limited. One of the few spectroscop ically confirmed 
cases is th e protocluster around 1 338 at z ~ 4. 1. IVenemans et all 
j2002l) and lVenemans et al, U2007I) have shown that the field around 
the radio gala xy is denser in LAE s than a blank field by a factor of 
4.8^ g. Also. lOverzier et al] d2006r) provided evidence for a rela- 
tivel y large number of Lym an Break Galaxies (LBGs) in this field. 

IVenemans et al] d2007l) found a total of 54 LAEs with EW > 
15 A in a field of 79.7 arcmin 2 around 1338. The narrowband fil- 
ter used was a custom filter with A c = 6199 A and a FWHM of 
59 A (hereafter NB620 for brevity), i.e. approximately twice as 
wide as the TF tunings used in this study. This difference in width 
implies that the same emission line, at the respective proper red- 
shifts, will yield a brighter magnitude in the TF. Applying the 50 
per cent magnitude limits to the 1338 catalogue will therefore not 
yield a proper comparison. The recovered Lycr flux is, however, 
relatively independent of filter width. We therefore use a cut of 
F Lya > 7.0 x 10~ 18 erg s" 1 cnr 2 . At z ~ 4.1 the width of NB620 im- 
plies a comoving volume of 12292 Mpc 3 . Applying the F Lya cut the 
density of LAEs i n the 1338 field is found to be 3.6 x 10~ 3 Mpc^ 3 . 

Applying the lVenemans et af, I J2007h selection criteria and ex- 
cluding the radio galaxy, we find a total of 5 unique LAEs with 
EWo > 15 A in the central protocluster field. With a volume of 
1570 Mpc 3 the density of candidate LAEs is thus (3.2 ± 1.4) x 
10~ 3 Mpc~ 3 where we used Poisson statistics for the lcr uncertainty. 
Thus the number density of LAEs around 6C0140 is comparable to 
that in the 1338 protocluster. 

The velocity distribution of the spectroscopically confirmed 
1338 LAEs have a FWHM of 625 ± 150 km s 1 . This is very narrow 
with respect to local galaxy clusters, but it is in agreement with 
the trend of decreasing velocity dispersion with increasing redshift 
dVenemans et alj|2007t) . 

The estimated redshifts of the candidate LAEs are compared 
to the expected redshift range of the protocluster in the lower panel 
of Fig. [6] Note that the location of the protocluster region does not 
coincide with the redshift of the radio galaxy. Instead we have cho- 
sen the redshift of the protocluster such that the number of LAEs 
consistent with being in the protocluster is maximised. This results 
in an offset with respect to the radio galaxy of Az = 0.0084 or 
Ai' ~ 465 km s _1 . This is consistent with what is observed in 
the 1338 protocluster where the 1338 radio galaxy is redshifted 
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Figure 6. Lower panel: the estimated redshifts of the individual objects with 
their respective uncertainties. The dark shaded region represents the redshift 
interval the protocluster would have based on the 1338 velocity distribution, 
whereas the light shaded region takes into account the l<x uncertainty. The 
protocluster region is chosen such as to maximise the number of protoclus- 
ter candidates. The spectroscopic redshift of the radio galaxy (z = 4.413) is 
marked by the vertical dashed line. The dotted line at z ~ 4.38 indicates the 
approximate redshift limit set by the edge of the circular protocluster field as 
shown in Fig. [3] Upper panel: distribution of the estimated redshifts. Also 
shown is a curve that indicates how the effective selection area varies as 
function of redshift. 

by 440 km s (Az = 0.0075) with respect to the majority of the 
confirmed line emitters. Thus the radio galaxy does not have to be 
at the centre of the structure in redshift space. In the situation as 
shown in Fig. [6] a total of 10 LAEs are consistent with being in 
the protocluster. This number decreases to 4 when we assume that 
the protocluster is centred on the radio galaxy. Note that a simi- 
lar displacement of the radio galaxy with respect to the bulk of the 
galaxies is seen for the spatial distribution of both the 6C0140 and 
1338 fields. Both radio galaxies are located not at the centre of the 
spatial distribution of emitters, but more at the edge. 

Also shown in the upper panel of Fig. |6]is the distribution of 
the estimated redshifts and a curve that indicates how the effective 
selection area varies as function of redshift. The area was estimated 
by determining the portion of the image for which A\ ow < A < Ahigh 
with A low and A hiai being respectively the lower edge of the TF1 
tuning and the upper edge of the TF3 tuning. The curve indicates 
that, based on the effective selection area, we would expect the ma- 
jority of the objects to have 4.35 < z < 4.39. However, we find 
a disproportionately large number of objects at z > 4.38 indicating 
that there is some concentration of LAEs close to the redshift of the 
radio galaxy. Applying a Kolmogorov-Smirnov test we determine 
that there is a probability of 4 x 10~ 3 that the observed distribu- 
tion is drawn from the expected distribution. The two distributions 
therefore differ at the ~ 3<x level. 

The top panel of Fig. [6] also shows that we are unable to ob- 
serve objects that are located at z > 4.42. This makes the reported 
overdensity of 8 ± 5 difficult to interpret. It may be that the proto- 
cluster structure extends beyond z > 4.42. If this is the case, then 
the true overdensity may differ from the value presented here. Like- 
wise, if the distribution of A c of the TFs had been chosen to probe 
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larger redshift values, then such a blueshifted overdensity as found 
for 6C0140 may be underestimated or even missed altogether. 

How do our results hold up when we account for the estimated 
contamination fraction discussed in Sect. 13.31 ? In both the compari- 
son with a blank field and the comparison with a z ~ 4.1 protoclus- 
ter we found 8 emitters in the 6C0140 field. Based on the minimum 
success rate of 75 per cent we therefore expect two interlopers in 
our 'protocluster' sample. Redoing the comparison with a sample 
of six emitters the followin g results are obtained . In the compar- 
ison with the blank field of iDawson et alj J2007I) it is found that 
the 6C0140 field is denser by a factor of 7 ± 4. Thus the 6C0140 
field harbours an overdensity of S g = 6 ± 4. The corrected number 
density in the 6C0140 field is (2.5 ± 1.3) x 1CT 3 Mpc~ 3 , which is 
also still in agreement with the 1338 field. The results presented 
here are therefore valid when contamination is taken into account 
and we conclude that the 6C0140 field is similarly overdense as the 
1338 protocluster. This indicates that it may evolve into a massive 
galaxy cluster at z = 0. Furthermore, this result supports the hy- 
pothesis that HzRGs are good tracers for galaxy overdensities in 
the early Universe. 



5 CONCLUSIONS & OUTLOOK 

We have presented the first search for high-z protoclusters employ- 
ing tunable narrowband filters. This pilot study focuses on the ra- 
dio galaxy 6C0140+326 at z ~ 4.4. Using a combination of three 
TF tunings we find a total of 27 unique LAEs in the field around 
6C0140+326. Division of the field in a protocluster and a fore- 
ground field shows that the protocluster field contains significantly 
more objects than the foreground field. This indicates that there is 
a concentration of LAEs near the redshift of the radio galaxy. 

A comparison to a blank field shows that the 6C0140 proto- 
cluster field contains an overdensity of 8 + 5. The number density in 
the protocluster field is also comparable to that found in the 1338 
protocluster at z ~ 4.1. Both these results are robust when taking 
into account the possible presence of interlopers. 

With the availability of three separate TF tunings we also es- 
timate the redshift distribution of the LAEs. Using results obtained 
for the 1338 protocluster we find that 4-10 of the LAEs have red- 
shifts consistent with being in a redshift interval spanned by a typ- 
ical z ~ 4 protocluster. Also, the redshift distribution is different at 
the 3cr level from the expected distribution with a relatively large 
number of objects at z > 4.38. This further strengthens the notion 
that there is a concentration of LAEs near the radio galaxy. 

These results are further evidence that HzRGs pinpoint high 
density regions in the early Universe. The overdensity around 
6C0140 may collapse at a later time to form a structure similar to a 
local galaxy group or cluster. Spectroscopic follow-up is needed to 
confirm this result. 

We have shown that tunable filters are an excellent method of 
confirming the presence of protoclusters around HzRGs at any red- 
shift. At the moment the wavelength range accessible to the red TF 
used in this study is limited to A > 6500 A and therefore z > 4.3. 
However, a blue TF covering the wavelength range A < 6500 A will 
be commissioned in the near future. This will open up the redshift 
range 2 < z < 4 which is where most of the known HzRGs are 
located. Our allocated GTC ESO large programme can then sig- 
nificantly expand the sample of protoclusters across cosmic time 
and this would allow an in-depth study of the evolution of these 
structures. 
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